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I. INTRODUCTION 

Over the past few years, there has been considerable 
interest in the study of small scale systems, with special 
emphasis on the issue of efficiency ;il . It was shown, 
on the basis of general thermodynamic arguments @ 
and using stochastic thermodynamics ||, that the effi- 
ciency at maximum power rf of a thermal engine, op- 
erating between a hot and cold bath at temperatures 
T^ 1 ) and T( 2 \ respectively, possesses universal proper- 
ties when expanded in terms of the Carnot efficiency 
?7C = 1-T (2) /T (1) : 
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This result is valid for "strong coupling," meaning that 
the particle and energy fluxes are proportional to each 
other. The value 1/8 for the coefficient of the quadratic 
term in addition requires a left /right symmetry, i.e., re- 
versal of fluxes upon inversion of forces. The above uni- 
versality ultimately derives from the reversibility of the 
underlying microscopic laws [Ijj. It has been verified 
in various models including transport of electrons 
through a quantum dot Q and of photons in a maser 
model Q . The purpose of this brief report is to present 
an analogous calculation for the transport of classical 
particles. As expected, universality is reproduced. An 
additional benefit is that, contrary to the case of quan- 
tum transport, an explicit analytic expression is obtained 
for 77*, namely, 
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II. MODEL AND MASTER EQUATION 

We consider a (small) reservoir of classical non- 
interacting particles which for simplicity all have the 
same energy e. This small reservoir constitutes our sys- 
tem. The number of particles in the system will be de- 
noted by n. This system exchanges particles with sev- 
eral other particle reservoirs v with temperatures T^ v > 
and chemical potentials p^ v \ respectively. We assume 



that the exchange can be described by Markovian dy- 
namics. This is the case, for example, if the transitions 
between system and reservoirs correspond to thermally 
activated processes over sufficiently high potential barri- 
ers. The probability distribution p n (t) for the system to 
be in state n at time t thus obeys the following master 
equation (n 6 N, quantities with negative n-index being 
zero by definition): 
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where W n +i^ n and W n —i, n are the rates (probabilities per 
unit time) for transitions n — > n+1 (the system gains one 
particle) and n — > n — 1 (loses one particle), respectively. 

From here on we are interested in the steady state op- 
eration of our system. The corresponding probability 
distribution p st is determined by the set of equations 
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or (since the flux is zero at the boundaries n 

n = 00), 
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= 0, 

(4) 
and 



(5) 



Note that the one-step hopping dynamics of our mas- 
ter equation has the peculiarity that the above formal 
condition of detailed balance is satisfied, at least with re- 
spect to the total transition rates W. This however does 
not necessarily correspond to true equilibrium, which re- 
quires detailed balance at the level of each of the separate 
processes taking place in the system, as we will see below. 

To proceed further, we have to specify the transition 
rates. We first assume that the exchanges of particles 
with different reservoirs are independent, so that the cor- 
responding rates W^ v ' add up to the total rate W, that 
is, 
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(6) 
(7) 



Second, statistical mechanics imposes physical con- 
straints on the separate rates Let us suppose that 
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the contact is broken with all reservoirs except for reser- 
voir v. The stationary state should in this case repro- 
duce the equilibrium state of the system in contact with 
this reservoir, p st — p e i^ v \ As is well know from equi- 
librium statistical mechanics, this is the so-called grand 
canonical distribution, which for classical non-interacting 
(ideal) particles is a Poisson distribution [gj], 
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The average particle occupation while in contact with 
reservoir v, fv- v \ is given by 



(9) 



where we have introduced the dimensionless quantity x Vl 

x v = /3l"\e- l iM), (10) 

and fjW = 1/ksT (ks = Boltzmann's constant). One 
can rewrite this equation in the more familiar form fj,^ — 

The requirement that p st = p e i^ for W = W {v) leads 
to the following "genuine" condition of detailed balance 
for these rates with respect to its equilibrium distribu- 
tion: 
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In the following we will adopt the standard choice of tran- 
sition rates Q that satisfy this condition (cf. law of mass 
action), namely, 
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with the n-independent rates obeying the balance condi- 
tion 



III. THERMAL ENGINE 
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Having identified the thermodynamically correct ex- 
pressions (fT2"j) and (ITUl) for the transition probabilities, 
we can proceed to a stochastic thermodynamic analysis 
(see for a brief review) of a system coupled to several 
reservoirs, with different temperatures and chemical po- 
tentials. One easily verifies that the steady state solution 
of ((3]), is again a Poisson distribution, 



Pn = rr e 



(14) 



with a steady state average number of particles reflecting 
the influence of each reservoir [compare with (|13[) ]: 
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This is most easily demonstrated by showing that the 
generating function is given by J2 n sn Pn = e^ -1 ^™ 

Even though Poissonian, this distribution corresponds 
to a nonequilibrium steady state. Using (TT2"|) -(fT5 j) . we 
find the following explicit expressions for the separate 
particle, energy, and heat fluxes from each reservoir v 
into the system: 
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Note that the fluxes from each reservoir to the system are 
strongly coupled, i.e., energy, heat and particle flux are 
proportional to each other. Furthermore, the fluxes from 
reservoir v are only zero when n = fv- v \ i.e. when the 
steady state distribution is the equilibrium distribution, 
cf. JSJ an d OS ■ The above Poissonian steady state (fT4|) 
does not, in general, obey detailed balance with respect 
to the separate rates W^ v \ and implies the following non- 
zero entropy production: 



s t = k B j2Y,( w > 



lPn- 



W^P^ln- 



W 



(") 



lPn 



5>J° 



Et(») y(v) 



k^fl] 



>o, 



In 



(19) 



where we have introduced the thermodynamic forces: 
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Of particular interest to us is the situation in which a 
heat current from a hot to a cold reservoir is used to drive 
particles uphill from low to high chemical potential. To 
investigate this case in more detail, we henceforth focus 
on the case of only two reservoirs v — 1,2, with reservoir 
1 the hot reservoir and 2 the cold one, > T {2 \ Note 
that we have not taken the thermal energy of the particles 
into account as this would require the consideration of 
a third thermal reservoir, making the comparison with 
Carnot efficiency more involved. At the steady state, one 
finds, using © and (fT3"|) . the following explicit results for 
the fluxes: 
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where we have introduced the rate 
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The entropy production reduces to the simple expres- 
sion 



Si = ksK(x2 — Xx)(e Xl 



> 0. 



(25) 



The power V of the engine, being the amount of net 
chemical energy produced per unit time, is given by 



v = ( M2 - mOjW 



(26) 
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The corresponding efficiency reads: 

V = —) — = 1 - I 1 - ^c) 
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(27) 



Before turning to the main issue of efficiency at maximum 
power, we first note that equilibrium, i.e. zero entropy 
production, cf. (|25p . is attained when x\ = x 2 . This 
does not require that the temperatures and T^and 
chemical potentials and iS 2 ^ be separately equal, a 
feature which is well-known for strongly coupled systems 
Q. In the vicinity of such a point, the machine can oper- 
ate reversibly, and its efficiency attains Carnot efficiency, 

T] = VC, Cf. (EH). 

Let us now turn our attention to the point of maximum 
power. From 
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one finds: 
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The corresponding efficiency reads: 



„* = at = a£ + *k + (32) 

1 V c-(l-Vc)Hl-r,c) 2 + 8 ^ •"' 1 ' 

which displays the expected universality announced ear- 
lier. 



IV. DISCUSSION 

We close with a number of additional comments. The 
above result is identical to the one obtained for a model 
based on particle transport via Kramers' escape [Io| . 
This can be understood from the fact that our model 
reduces to this case when the outgoing rates fc_ become 
very large. In this limit the number of particles in the 
system goes to zero, and the only remaining processes 
are the thermally activated transitions from one reservoir 



into another, via fast passage through the system, which 
plays the role of a short-lived activated state. The fact 
that our more general model reproduces the same result 
as Kramers' escape suggests that ^ may have a wider 
applicability in classical transport. The above results for 
particle flux, power, entropy production and efficiency at 
maximum power are also reproduced by taking the clas- 
sical limit in the problem of electron transport through 
a quantum dot (5). 

The above model for particle transport can be repre- 
sented as a simple chemical reaction, namely, 



(33) 



where n represents the number of particles of species 
X, which can transmute into the species JfM whose 
chemical potential is fixed by reservoir constraints. This 
representation is most natural in an isothermal system, 

Finally, we mention a similarity of the model studied 
here with that of an underdamped Brownian particle in 
contact with several heat baths. The Langevin equation 
for such a particle has the following form: 



mv 



(34) 



where are independent normal white noises. The 
noise intensity is chosen in accordance with the 
fluctuation-dissipation theorem so that the stationary 
distribution when in contact with each reservoir sepa- 
rately reduces to the corresponding Maxwellian veloc- 
ity distribution. As in the model studied here, the sta- 
tionary distribution of the system in simultaneous con- 
tact with multiple reservoirs has an " equilibrium shape" , 
i.e., a Maxwellian distribution, however at a temperature 
which is the geometric mean of the bath temperatures. 
Detailed balance is broken at the level of the exchange be- 
tween particle and separate reservoirs, with the Brownian 
particle functioning as a thermal contact between reser- 
voirs. This Brownian model has been studied in great 
detail, revealing detailed properties of the corresponding 
nonequilibrium steady state [Tl| . The model presented 
here has the additional advantage of allowing both heat 
and particle transport. We therefore expect it to be an 
interesting candidate for revealing further properties of 
nonequilibrium steady states. 
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